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Abstract 

In this paper we study the Einstein-Boltzmann system with Bianchi I 
symmetry. We show that for small initial data the corresponding solutions 
of the Einstein-Boltzmann system are future geodesically complete and 
that they isotropize and have a dust-like behaviour at late times. Detailed 
information about the metric and the matter terms is obtained, and the 
results show that the solutions tend asymptotically to the Einstein-de 
Sitter solution. 


1 Introduction 

In theoretical cosmology one wants to understand the future and the past be¬ 
haviour of different cosmological models. This is in general too ambitious: in¬ 
stead, one looks at simplified models, such as the important class of homoge¬ 
neous spacetimes. Before obtaining a global and detailed picture of the future 
and past behaviour of a cosmological model, one hopes to obtain an under¬ 
standing of the asymptotic behaviour, i.e., what happened close to the initial 
singularity (the big bang) and what is the distant future of our Universe. Both 
problems are still far from being trivial. Indeed, only very special models are 
understood to date, and the predictions are very different from each other. In 
this paper we will only deal with the future direction. 

How do the asymptotics depend on the way the matter content of the Uni¬ 
verse is described? It is standard to model the Universe with a perfect fluid, 
but this is only one of the several options for the choice of the matter model 
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that exist, and it is of interest to know if and how the dynamics change with a 
different matter model. The most common ones apart from the perfect fluid are 
elastic matter, field-theoretic models and kinetic models. We are interested in 
kinetic models and will focus in the present paper on the case where no cosmo¬ 
logical constant is present. It is mathematically more challenging and we plan 
to investigate the cosmological constant case in the future. In the last years 
there has been a lot of progress concerning the dynamics in the Einstein-Vlasov 
case. For an overview of different results we refer to [^. There the dynam¬ 
ics of collisionless matter is analyzed for Bianchi spacetimes which are locally 
rotationally symmetric (LRS). In absence of LRS symmetry, the asymptotic 
behaviour in the case of homogeneous spacetimes has been studied for a large 
class of solutions (cf. [30). On the other hand, not so much has been done 
in the Einstein-Boltzmann case. One of the reasons is that the formulation of 
the Einstein-Boltzmann system is not trivial in general relativity (cf. [i,0). 
Instead, one often considers certain symmetries on the spacetime to s implify 
the full system of equations, and can find several results [l|, [3[3l3,l3lii 
for Bianchi type I case. In the present paper we also consider the Bianchi type 
I spacetime, and the objective of this paper is twofold. First, we study non¬ 
diagonal Bianchi I spacetime. In the results mentioned above the authors only 
deal with the diagonal case, for instance FLRW or LRS cases, but our result is 
valid for spacetimes of Bianchi type I which are not necessarily diagonal. Sec¬ 
ond, we obtain asymptotic behaviour of solutions at late times. For instance in 
0 , the authors considered large initial data and only obtained global existence 
of solutions, but in this paper we consider small solutions and obtain late time 
behaviour. 

This paper is organized as follows. In Section [2l we derive the Einstein- 
Boltzmann system with Bianchi symmetry. We first review the Einstein-Vlasov 
system with Bianchi symmetry, and then parametrize the mass shell with an 
orthonormal frame to obtain the Einstein-Boltzmann system with Bianchi sym¬ 
metry. In Section |3] we present the equations we will work with, namely the 
Einstein-Boltzmann system with Bianchi I symmetry. In Section SI we collect 
several elementary lemmas, which are used to estimate the collision operator of 
the Boltzmann equation. Global-in-time existence of solutions to the Einstein- 
Boltzmann system will be proved by a standard iteration method. In Section 
SI we assume that a distribution function is given, and then show that clas¬ 
sical solutions to the Einstein equations exist globally in time. In Section SI 
we assume that a metric is given, and then show that classical solutions of the 
Boltzmann equation exist globally in time. We combine the results of the pre¬ 
vious sections to obtain the main results. Future geodesic completeness and 
the asymptotic behaviour of solutions to the Einstein-Boltzmann system with 
Bianchi I symmetry are obtained under suitable small data assumptions. De¬ 
tails on the relativistic Boltzmann equation in the Minkowski case are collected 
in the appendix. 
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2 The Einstein-Boltzmann system with Bianchi 
symmetry 

In this section we derive the Einstein-Boltzmann system with Bianchi symmetry. 
We introduce the Einstein-Vlasov system in the first subsection. In the second 
subsection we consider this system with Bianchi symmetry. Then parametrizing 
the mass shell with an orthonormal frame we obtain the Boltzmann equation 
with Bianchi symmetry in the third subsection. In the last subsection we present 
two different parametrizations of the post-collision momenta which will be nec¬ 
essary in the following. 

2.1 The Einstein-Vlasov system 

A cosmological model is described via a time-orientable Lorentzian metric 
(we will use signature (—,+,+,+)) on a manifold At, which means that at each 
point of At the two halves of the light cone can be labelled past and future in a 
way which varies continuously from point to point. This enables to distinguish 
between future-pointing and past-pointing timelike vectors. Using units such 
that the gravitational constant and the speed of light in vacuum are equal to 
one, the interaction between the geometry and the matter is described by the 
Einstein field equations: 


= SttT, 


a/3 j 


where Gag is the Einstein tensor and Tag is the energy-momentum tensor. 
Indices are lowered and raised by the spacetime metric gag and its inverse 
We use Greek indices for indices which run from 0 to 3 and Latin indices from 
1 to 3. The energy-momentum tensor describes the matter content, and we will 
use a kinetic theory to describe it. We refer to [i, 22| for an introduction to the 
Einstein-Vlasov system. 

Consider now a particle which moves only under the influence of the grav¬ 
itational field. The worldline x°‘ of a particle is a null or timelike geodesic 
depending on whether it has zero or non-zero rest mass. The tangent vector to 
this curve is the 4-velocity and pa is the 4-momentum of the particle. Let 
TxM be the tangent space and T*M. the cotangent space at a point in the 
spacetime M. For particles of the same type and with the same (constant) rest 
mass m, the mass shell relation is given by: 


PaPgg^^^ +rn^ = 0, 


where the Einstein summation convention that repeated indices are to be summed 
over is used. The component po can be expressed in terms of the other compo¬ 
nents via this relation: 

Po = ^ ( - ± \/{Pag°°-Y - g°°{PaPbg‘^^ + ■ 

We are interested in future pointing momenta, i.e., 

P° = Pog°° + Pag^"" = +\J (Pa5°“)^ - 5°°(PaP&g“'' + -m?) > 0. 
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Define the phase space Pm for particles of mass m: 


Pm = {(a;“,Pa) : x°' € M, Pa € TfM, PaPpg^'^ + = 0, p° > 0}, 


which is a subset of the cotangent bundle T*M = {{x°‘,pa) ■ x^ G M, Pa G 
TfM}. Usually the geodesic equations are expressed in terms of u“, but one 
can do the same with Va- Using the convention that V EpPa = ^fpE-y, V being 
the covariant derivative and Tapp = gaS^pp, we have that 


0 = V„ 


Xv^ep) = = (^-v^v^Tpsa + 


The geodesic equations are then: 




where t is an affine parameter. The components of the metric connection Tapp 
can be expressed via the Koszul formula in the vector basis Ea with commuta¬ 
tors [Ea,Ep] = 'y'^i^Ep, where 7 ^^ are called commutation functions: 



(ppig CXfl ) + Ep{gpa) 


Eaigpp) + ippgas + lipgps 



( 1 ) 


Using the fact that we are summing over v°‘v^ the geodesic equation for vp can 
thus be expressed as 


^ ^ ^ [Epigap) + 2jf,pgas] 

The restriction of the Liouville operator to the mass shell is defined as: 

^ dx°‘ d dpb d 

ds dx°‘ ds dpb 


Using the geodesic equations it has the following form 


L = P 




■ mTabpV°v'^ 


dpb 


This operator is sometimes also called geodesic spray. The matter is described 
by a nonnegative real valued distribution function f{x°‘,pa) on Pm- In the 
Vlasov case we have that 


L{f) = 0. 

It remains to define the energy-momentum tensor Tap in terms of the distribu¬ 
tion and the metric. Before that we need a Lorentz invariant volume element on 
the mass shell. The cotangent space has the Lorentz invariant volume element 
[— det 2 dpodpidp 2 dp 3 , where det is the determinant of the spacetime 
metric. Now considering po as a dependent variable the induced Riemannian 
volume of the mass shell considered as a hypersurface in the cotangent space at 
that point is 

w= detg''‘^^]~^dp^,, 

pv 
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where dp^ = dpidp 2 dp 3 . Now we define the energy momentum tensor as follows: 



One can show that Tap is divergence-free and thus it is compatible with the 
Einstein equations. For kinetic matter all the energy conditions hold. In partic¬ 
ular the dominant energy condition is equivalent to the statement that in any 
orthonormal basis the energy density dominates the other components of Tap, 
i.e., Tap < Too for each a, /3. Using the mass shell relation one can see that this 
holds for kinetic matter. The nonnegative sum pressures condition is in our case 
equivalent to gabT°‘^ > 0. Another important quantity is the particle current 
density, which is defined as 

iv“=y 

It holds that 

VaN^ = EaN^ + T^aN^ = 0 , 

which is an expression of the conservation of the number of particles. For 
simplicity we will assume that our system of particles is made of a single species 
of particles (no mixtures) and that all particles have unit mass, i.e., the units are 
chosen such that m = 1, so that a distinction between velocities and momenta 
is not necessary and the mass shell relation is 

PaPpg"^^ = -m? = - 1 . ( 2 ) 


2.2 The 3+1 decomposition of the Einstein-Vlasov system 
with Bianchi symmetry 

A Bianchi spacetime is defined to be a spatially homogeneous spacetime whose 
isometry group possesses a three-dimensional subgroup that acts simply tran¬ 
sitively on the spacelike orbits. A Bianchi spacetime admits a Lie algebra of 
Killing vector fields. These vector fields are tangent to the group orbits, which 
are the surfaces of homogeneity. Using a left-invariant frame, the metric in¬ 
duced on the spacelike hypersurfaces depends only on the time variable. If 
denote the one-forms dual to the spatial frame vectors, the metric of a Bianchi 
spacetime in the left-invariant frame is written as 

gU) = g = gab{t)uj°-{t)uj’^{t). (3) 


Due to the simple form of the metric in the Bianchi case, it is useful to do 
the 3-1-1 decomposition of the Einstein equations in the left-invariant frame. 
Using the sign conventions as in 21| we obtain with zero shift and the lapse 
function equal to one (/3 = 0 and a = 1 in (2.29) and (2.34) of nil) the following 
evolution equations: 


gab — ‘^^ah, 

^ab — Eab ^ ^ab ^(ll ^bd^^ac ^'^Sab ^'^gabi.S /?), 
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where kab is the second fundamental form, k = kabg‘^^ its trace, Rab the Ricci 
tensor of the spatial metric gab, S = g°‘^Sab, and the dot denotes the derivative 
with respect to time t. The constraint equations are given as ((2.26) and (2.27) 
of [HI) 

R — k°'^kab + k‘^ = 167rp, 

V“fcob = -STrTof,, 

where R = Rabg^"^ is the scalar curvature. Since we use a left-invariant frame the 
distribution function / will not depend on x°'. Moreover since goo = g^^ = 
and = 0 we have that = —po = Defining p = Too, 

Ja = —Toa and Sij = Ty the frame components of the energy-momentum 
tensor and the particle current density are thus 


P = (detg)-5 

/ f{t,P*)'J^+PaPbg°-^dp^: 



Ji = (detg)“^ 

/ f{t,p*)pidp^, 


JS.3 

Sij = (detg)“^ 

f fit,P*) n -^- -^dp: 

= (detg)-5 

/ f{t,P*)dp*, 


7r3 

V* = (detg)-^ 

/ fit,P*) n - -rdp: 

^R^ V 1 + PaPbg°-^ 


The conservation of the number of particles implies then 

iv°-b = 0. 

We now consider the Vlasov equation. Since b is spacelike 

'^c,b^l.P°‘P^ = llbdaSP^^P^- 
Consider now two different cases. 


(4) 


• In the metric approach, a basis of the frame vectors is chosen, such that 
they are time-independent and the commutation functions are made 
equal to the structure constants of the Lie algebra which are denoted by 
Cj^. Only the spatial components remain and as a result we have 


dt 


+ ^^abP‘'Pd 


K 

dpb 


= 0 . 


• In the orthonormal frame approach, gab = Pab- On the other hand the 
commutation functions and the spacelike frame vectors depend on time. 
However 7 °^ = 0 and 7 qj = k^, where we use a hat to indicate that we 
are in the orthonormal frame. As a result 

f+ (i?W + ()i")-Siw!i“)^=0. 


For more details for these approaches we refer to [l^- Now by coupling the 
Vlasov equation to the Einstein equations through the energy-momentum ten¬ 
sor, we obtain the Einstein-Vlasov system with Bianchi symmetry. 
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2.3 The Boltzmann equation with Bianchi symmetry 


We now have to introduce the collision operator. The idea is to use first the 
orthonormal frame, because then the collision operator has the same form as 
for the Minkowski case. Then, we have 


% + 


ifr^itbPdP 




where Q{f, /) is the collision operator as in the Minkowski case. If we make the 
transformations from one frame to the other we arrive at 

% + ^C^abP‘'Pd^=Q(.fJ), 


dt po 


dpb 


where Q{f, /) is now the collision operator for the Bianchi case and / depends on 
p^,. The details for the collision operator in the Minkowski case are presented 
in the appendix. Let us introduce now some standard terminology. If two 
particles with momenta p“ and collide, due to the conservation of energy 
and momentum the following relation holds between the momenta (p“, g“) and 
the post-collision momenta {p'°‘,q'°‘): 

-p q'°‘ =p°‘+ q°‘. 

The post-collision momenta will be parametrized in two ways with param¬ 
eters The relative momentum h and the total energy s are defined 

by 


h= ViPa- 9a)(p“ - g“), S = -(Pa + ga)(p“ + g“)- 
Using these variables the Mpller velocity vm is defined as 

hy/s 


VM = 


,0„0 ■ 


4p0qt 


For notational simplicity we also use 


n“ = p“ -H 9° 


u = p 


and then we may also write = UaU°‘ and s = —nanP. Finally we denote by 
a = a{h, 9) the scattering kernel, where 9 is the scattering angle, which is the 
angle between and u'°‘. For background on the relativistic kinetic theory and 
intuition on the variables and expressions introduced, we refer to and pages 
94-116 of For a more general presentation of the non-relativistic Boltzmann 
equation we refer to 261. 


2.4 


Two representations of the collision operator and the 
post-collision momenta 


We generalize the Boltzmann collision operator in the Minkowski case to the 
Bianchi case. One way to do this is to use an orthonormal frame in the momen¬ 
tum space 13|. Applying the usual Gram-Schmidt process, which also works 
for pseudo-Riemannian metrics as mentioned in [l^ , we obtain an orthonormal 
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frame {e^} from a given frame {Ea} such that = e^Ea or Ea = where 

is the inverse of e“. It satisfies the following relations: 


P“ = e“p'", gap = 9i^9'pr]f,„, 

where Pa = PafiP^ and P/^ = r/^^pE. In the Bianchi case with the metric ([3]) , we 

J. 
dt 


may take Cq = ^ to obtain 


e° = 0° = l, eS = e° = 0g = 0° = O, a =1,2,3. 
Thus, po = po, = P^, and 

^ nb ^ _ ncnd^ ab ^a^bcd 

Pa — ^aPb-) P — ^bP : 9ah — ^a^bPcd: 9 — ^c^dV 

We can also compute the other components as follows: 



66 = 


and its inverse 


/ 

y'ml 
—.912 


0 

.911 


.912.923 —.913.922 


v'iriTv'mi' 

-.911.923+912913 911922 —(912) 

\/m2\/ms A/moA/ms 



/ 911 

y/mi 

912 
V^rni" 
L 913 
\ ^mi 


0 


£12_ 911922 —(912) 

ym\ ^mi ^1712 

.913 911.923—912913 


0 

0 

det g 

y/m^y/rn^ 


(5) 


( 6 ) 


where m-i, m 2 , and m 3 are the principal minors of the matrix gab, i.e., the 
quantities gu, 311(722 — ( 312 )^) and det 3 , respectively. In fact, the matrix 9^ 
is uniquely determined, since gab is symmetric and positive definite, which is 
known as the Cholesky factorization. Note that in order to estimate the quan¬ 
tities and 9^ we need to estimate the lower bounds of the principal minors of 
the matrix gab- The following two lemmas give the lower bounds of the principal 
minors. Below, denotes the set of all n x n matrices over R. For the proof 
of the lemmas we refer to Q. 

Lemma 1. (Hadamard’s inequality) Let A = (aij) G Mn be positive definite. 
Then, 

det id dll ■ • * a^iji 

with equality if and only if A is diagonal. 

Lemma 2. (Fiseher’s inequality) Suppose that the partitioned symmetric matrix 


H = 


A B 
C 


G Mp^q, A G Mp, C G Mq 


is positive definite. Then, 


deti? < (deti4)(det C). 


The following estimates are now easily obtained, 

1 ^ 322333 f ^ 333 
mi “ det 3 ’ m 2 ~ det 3 ’ 


( 7 ) 














and we can conclude that the quantities and 9^ are bounded as long as 
(detg)“^ and gab are bounded. The above inequalities are also used to estimate 
the asymptotic behaviour of the components and 9^. 

In an orthonormal frame the representations of collision operator and post¬ 
collision momentum become identical to the Minkowski case (see Appendix). 
Hence, we have 


QifJ) = 


h^/s 


’{h,9)(j{p)f{q) - f{p)f{q)jdujdq. 


/b3 

Here, h and s are understood as 

h = - q>^)(p'' - q''), s =-r]f,^{p>^ + q>^){p''+ q"), 

where p^ and q^ denote the components of momenta in the orthonormal frame. 
The post-collision momentum is written by 


5/0 

~lk 


fP -I- (f 


pP ^q*" h 




h{n ■ Lo) 

2 yi 

fP 


(h ■ uj)r 


\ 


J 


where = p^ + q^. In Glassey-Strauss’s framework, the above quantities are 
written as follows: 


QifJ) = 

and 


hs{'hP)‘^a{h, 9) 


/s2 4p0g0((h0)2 - (h 


g)2)3/2 {fip')fiq') - fip)f{q))d^dq, 


g/O 


f p° + q° ^ h (h-g) \ 


V 


p'" + q'^ h 


2 ^(h 0 ) 2 -(h- 0 ^ 

2 y^(h 0 ) 2 - (h^ / 


We now switch to the p* and g* variables through the transformations e“ or 
9^, and then obtain representations of the collision operator and post-collision 
momentum in the original frame. In the end what we obtain is the collision 
operator Q in Strain’s framework, 

Q(/,/) = (detp)"5 / [ VM<j{h,9)(f{p'^)f{ql)-f{p^)f{q^))dujdq^ 

Jr3 y §2 V / 

with the parametrization of post-collision momenta 

/ P° + q\ hnaeluj^ \ 


p'O 

P^ 


V 


Pi + qi 


2 yi 
0 


2 2 

and in Glassey-Strauss’s framework 


+ o gia^b^ + I —p — 1 




g^’^ncUd 


QifJ) = (detp) 2 f [ 
JR3 Js- 


§2 ((n0)2 - (jiae 


3/2 (/(P*)/(g*) - f{P*)f{q*))didq, 
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with the parametrization of post-collision momenta 



/ / + 9° h ngete \ 

1 

o 

2 2 y'{n^Y - (nce55"*)2 

\ P'i ) 

Pi + di h n^giael^'^ 


\ 2 2 y^(n0)2 - (nce(jC'^)2 / 


For simplicity we used upper indices for zeroth components since they only 
differ in sign, i.e., Pq = —p'^, and moreover we will regard as well as as 
a function of spatial components through the mass shell condition ([2]) as in the 
Vlasov case. Note that the collision operator can be written as Q = Q+ — Q- 
in both cases, which are called the gain and the loss terms, respectively. 


3 The Einstein-Boltzmann system with Bianchi 
I symmetry 


The Einstein-Boltzmann system with Bianchi I symmetry is obtained by com¬ 
bining the equations derived so far. At the end of the section we make assump¬ 
tions on the scattering kernel and the Hubble variable. Several basic estimates 
of the Hubble variable will also be given. Note that the Ricci tensor, the scalar 
curvature and the coefficients of the metric connection with spacelike indices 
vanish in the Bianchi I case. Now, the evolution equations of the Einstein 
equations are 


Qah — 2 /^ 05 , (8) 

kab = {9‘'‘^kcd)kab-‘2{g‘'‘^kM)kac-8TTSab+^T^9ab{S - p), (9) 


and the constraint equations are 

- k°'^kab + = 167rp, 

0 = -SirTob. 


( 10 ) 

( 11 ) 


The first expression of the Boltzmann equation is 

^ = {det 9 )~i J ^ J^^VMcr{h, 0 )(^f{p'^)f{q[) - f{p^)f{q^)^dujdq^, (12) 

where 


/o 


P 

Pi 


Pi + 9i 


V 


+ 2 1 Sia^bU} + 


h rtgega;^ 

2 

VS 


ngCbOJ Uj 

9‘^‘^ncnd 


(13) 


and the second expression is 


dt 


= idetg)-i [ j 
jR3 JS 


VM\/s{n°)'^a{h,6) 


S 2 ((n 0)2 - (nae“C*’)^ 


j^(/(p'J/(9() - f{p*)f{q*))d^dq^, (14) 


where 



/ p° + q° h ngel^^ \ 

1 

O 

2 2 y'(n0)2 - (nce(j^'')2 

V p'l ) 

Pi + qi h 'nPgigel^^ 


\ 2 2 y'(n0)2 - (n,eS^<^)2 / 


( 15 ) 
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( 16 ) 


They are coupled to each other through the energy-momentum tensor 
P = (detg)"W /(t,p*)(l 

Sab = (detg)“W fit,p^)paPbii + g‘"^PcPd)~^dp^. (17) 

JR^ 


In this paper the Einstein-Boltzmann system with Bianchi I symmetry will refer 
to the system of equations ®-(ini), and global-in-time existence and asymptotic 
behaviour of solutions will be studied. We expect the spacetime to behave as 
the Einstein-de Sitter model at late times. Hence, it will be useful to use the 
following notation throughout the text for a scaled version of the spatial metric 
and its inverse: 


gab = t*gab and g°‘'" = t 


(18) 


Throughout the paper, and will denote the components of an orthonormal 
frame and its inverse, respectively. Change of variables between p and p will be 
used frequently. Indices will be lowered by gab for p“, i.e., Pa = gabp’’, but will 
be lowered by pab in the orthonormal frame case, i.e., Pa = PabP^■ 

In the rest of the section we make two assumptions on the Einstein-Boltzmann 
system. The first one is on the scattering kernel of the Boltzmann equation, and 
the second one is on the Hubble variable which will be defined below. Several 
basic estimates of the Hubble variable will also be given, and the estimates hold 
for any classical solutions. Those estimates will be improved in later sections 
by assuming smallness of initial data. 


3.1 Assumption on the scattering kernel 

The scattering kernel is a quantity which depends on the relative momentum h 
and the scattering angle 9. The relative momentum h is defined as = UaU°‘, 
and the scattering angle cos 9 = Uau'°‘/h‘^, where we used = UaU°‘ = 

Eor a technical reason, in the present paper, we will assume that the scattering 
kernel is given by cr = cr(p*, q*,u}). To be precise, we assume that the scattering 
kernel, when written in an orthonormal frame, i.e., 

CT = d{p,q,uj), 

satisfies the following: 

• It is non-negative and bounded. 

• It is differentiable with respect to p, and its derivatives are bounded. 

• It is supported in ||up < Irt'P -I- for some positive number C. 

The third assumption on the scattering kernel may not seem to be physically 
well-motivated. However, we remark that for any covariant variables p* and 
< 7 * the support will eventually cover all a; G By the estimate (|28ll we have 
|up — |u'p < h?{p^q^ — 1) and note that 

p° = V^ +g'^’^PaPb, = -(p° - q°f + g°-^{pa - qa){pb - qb)- 

7 4 

In the present paper we will show that g ^ t“3. This implies that the non¬ 
negative quantity h?(jfiq^ — I) decays to zero with a rate t~3. Hence, we can 
see that the inequality |up < \u'\'^ + Ct~3 will eventually hold. 
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3.2 Assumption on the Hubble variable and time origin 
choice 


Instead of the mean curvature, i.e., the trace of the second fundamental form 
k = g^-^kab, in cosmology often the Hubble variable is used: 



In this paper we assume that the Hubble variable H is initially positive (the 
trace of the second fundamental form is initially negative), i.e., 


H{to) > 0, 


which corresponds to an initially expanding universe. Then, the first constraint 
equation implies that if the spacetime is neither Bianchi IX nor flat, then H > 0 
for all time. The initial data set corresponds to the triple ,/°). Since 

one can always add an arbitrary constant to the time origin, we choose our time 
coordinate and time origin such that to = —2(fc°)“^, which is just some positive 
number. In other words we can set without loss of generality 

H{to) = (19) 


3.3 Basic estimates of the Hubble variable 


It is useful to express the second fundamental form in terms of its trace free 
part (Tab and the Hubble variable: 


^ab — (^ab Hgab- 

The evolution equation for the Hubble variable is 

• 47r 

H = -3H^ + —{3p- S). 

Let us use the notation 


F = 


1 

4IP 




( 20 ) 


( 21 ) 


to denote the shear. If we substitute the energy density with the constraint 
equation, then we obtain F = ^{1 — fl) with H The evolution equation 

for H can then be expressed as 


d{H-^) 3 ^ 47rS' 

—i- = — l- F -\ - 

dt 2 3iJ2 


( 22 ) 


From m using the inequality p> S and from (1^ using the fact that both F 
and S are positive we obtain: 



d{H-^) 

dt 


<3, 


which integrating leads to 



< H < -t 


-1 


(23) 
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where we used ©1). We remark that the estimate (©1) holds for any classical 
solution, and this estimate will be improved by assuming smallness. 

We collect several evolution equations which will be used to obtain detailed 
asymptotic behaviour in later sections. The evolution equation for F is 




47r 



(24) 


For the Bianchi I case the conservation of the number of particles (|3]) is simply 

N° = -3HN°. (25) 

Another useful relation concerns the determinant of the metric g ((2.30) of [llj): 


d(log detg) 
dt 


= -2k = 6H. 


(26) 


4 Estimates of the collision kernel 

In this section we collect several elementary lemmas concerning estimates of 
quantities related to the Boltzmann collision operator. The proofs are straight¬ 
forward, and the lemmas will be used in later sections to estimate the distribu¬ 
tion function of the Boltzmann equation. 

Lemma 3. The following estimates hold: 


s = 4 -|- h"^, 2 < y/s, 

h < y/s < 2y/p^q^, 

(27) 

h? < UaU°' < p^(fh^, 


(28) 

9° Q.n 

\ ^ < s < (n° 2 - 

pU grU 

Me?. 

(29) 


Proof. It is easy to compute that 

s = A + h'^, 

which obviously gives us that 

^/s > h, -v/s > 2. 

For the inequalities in (l?7l) and (1^ . we use 

(1 ± gabP°‘q^Y = 1 ± '2gabP°‘q^ + {gabP°'q^)'^ 

< 1 + + gabq^q^ + {gabpy^){gcdq^q'^) = {p°)Hq°?- 

Then, the last inequality of (EZD is verified as follows: 

s = -{pa + qa){p°' + q°‘) = 2 + 2p°q° - 2gabP°’q’' < Ap°q°. 
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Since h? = — + UaU°‘, the first inequality of (E51) is clear, and the second 

inequality is proved in a similar way: 




= 2p°q° - 2(1 + gijp'q^) 

' V p°q° + 1 + QcdP'^q'^ 

1 + gabP>'’ + 9abqy^ + {9abP>'’){9cdq^q’') - 1 - ^g^bP^q^ - {gabP'^q^? 


pOqO 


pOqO 


> 


gabip^^ - q^W - q^) 


pOqO 


where we used the fact that gob = 0 . 

For the estimate (1^^ . we use the definition of s as follows: 


s = {n°f - gab-n^'n!’ 

= {n^f - {gabn^n^){gcde‘rCe'^^n < - {gabU^elO^ 

where we used gabS^e^ = Pcd, and this gives the upper bound of s. The lower 
bound for s in ( 1 ^ is proved as 


s = 2 + 2 p°g° - 2gabp'"q^ > 2 + 2p°q° - 2\/gabP°^P^\/gcdq^'q^ 


o I J ip°)Hq°? - 9abPVgcdq‘^q^ \ 

V p°g° + \/gabP°-p^\/gcdq’^q’^ J 


>2 + 


p°g° + \/gabP°-p^\/gcdq'^q'^ 

1 + 5abP>'' + gabq’^q^ 2p°q° + {p°)^ + {q° f - 1 ^ g° p° 

— „0 „0 ’ 


pOqO 


pOqO 


and this completes the proof of the lemma. 
Lemma 4. The following estimates hold: 


pu qu 


□ 


a 0 P 

dp,P =^, 

h [p° gO 

^ q° (p^ g* 

dp, Vs = ^ 

ys Vp q 

dpi\J{n°Y - (naelff>f = 


yJ{n°Y - (uaeli^y 


^.,.(P^ qgeteelV 

19 qO 


,ofV Pae^b^elV 
^ \p° pO 


dr, 


naelw^nj 


Jlrn^ 


el^uj°'nj ^ nae1,u!^6j ^noe^a;*’njg®'^nc 
; 2 ,- 


711)71^ 


TlrTl^ 


{ndU’^Y 


(30) 

(31) 

(32) 

(33) 

(34) 


Proof. The results of this lemma are direct calculations. For instance, dSOD is 
obtained as follows: 


9piP° = 9p. v'l +PaP“ = 






vT+p^ p° ’ 
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where <5^ is the Kronecker delta. The other identities are similarly obtained, 
and we skip the proof. □ 


Lemma 5. Let p* and q^, be given. Suppose that and q'^ are post-collision 
momenta with a parameter cu which satisfies the assumption on the scattering 
kernel. Then, we have 

9°'^(,PaPb + QaQb - p'aP'b “ < C'. 

where is defined by m- 

Proof. We first consider g°‘^ to take an orthonormal frame such that 

5“'’(PaP6 + qaqb - p'aP'b “ = \P? + “ \p'? “ W ? ■ 

Then, we use the energy conservation of particles to have 

\p? + \q? - \p? - \q? = {py + y? - {p’y - yy = y°q° - 2 / 9 °. 

Since h is an invariant quantity in the collision process, we have 

= (Pa - 9a)(p“ - 9 “) = (Pa “ O (p'“ “ 9'“), 

and then 

h^ + y + qy = \y+ 

By the assumption on the scattering kernel, we have 

4p'°g'° - 4/g° = \y - < Ct-i. 

Since g°‘^ = t~ig°‘^, this concludes the lemma. □ 

5 Existence of solutions of the Einstein equa¬ 
tions 

In this section we study the Einstein equations for prescribed matter terms. We 
will assume that a distribution function / is given, and the energy-momentum 
tensor is derived from (HU) and (O, and will show that classical solutions to the 
Einstein equations (©-dHl) exist globally in time and have certain asymptotic 
behaviours. In this section we assume that the distribution function satisfies in 
an orthonormal frame the following estimate: 

f{t,p*) = f{t,p) < C'/exp(-t3|pp), (35) 

where (7/ is a positive constant which will be determined later Below, and in 
the rest of the paper, C will denote a universal constant which depends only on 
initial data and may change from line to line, but the constants with subscripts, 
such as Cf, will be understood as fixed ones. 


15 


5.1 Global existence of solutions 

Global-in-time existence of solutions to the Einstein equations with a given 
distribution function is easily proved by following the arguments of [^ . where 
the Einstein-Vlasov system has been studied for several different Bianchi type 
symmetries. In the Vlasov case, it is usual to assume that the distribution 
function / has a compact support, i.e., there exists a positive number P such 
that if \pa\ > P^ then f(t,Pf) = 0. In this case, we can estimate the matter 
terms and (113 as follows: 

P < AP\ < AP\ 

where A is a positive constant depending on (det^)”^ and pab- Hence, as long 
as (detg)”^ and pab are bounded, the matter terms are also bounded. Itera¬ 
tion scheme then works well, and local-in-time existence follows. For detailed 
arguments we refer to (20l| . In the Boltzmann case, we do not assume that the 
distribution function has a compact support, but instead assume that / decays 
at infinity. If we assume dsni), then the matter terms are estimated as follows: 

P= ( f{t,p){^ + \p?)^d'P<Cf [ exp(-<3|pp)(l-P |pp)5dp < 

JR3 

where p is written in an orthonormal frame. On the other hand, we use Pa = 
9 abe^p‘^ = Vcbd^P'^ to estimate Sab as follows: 

Sab = {detp)~i f{t,p*)paPb{l+g‘'‘^PcPd)~^dp^= / f(t,p)paPbil + \p\‘^)~^dp 
jR3 JR3 

= [ f{t,p)pcPd{^ + \p\'^)~^dp = e^e^Scd, 

JS.3 

where Pc = PacP°‘, and Scd is easily estimated: 

|-5cd| < [ f{t,P)\Pc\\Pd\{l+\p\'^r^dp<Cf [ exp(-t“|pp)|pp# < 

JR3 JR3 

_ (36) 

By the estimates (0, the quantities 9^ are also bounded as long as (det p) ^ and 
Pab are bounded, and we conclude that the matter terms are bounded. Applying 
the same iteration as in [^ . we obtain local-in-time existence of solutions pab 
and kab to the Einstein equations. Now, we can extend the local existence to 
global existence by the arguments exactly as in [^, and obtain the following 
global-in-time existence: 

Lemma 6. Suppose that f = fit,p) is piven such that it is and satisfies 
(IMll . If Pab{to) and kabitfi) are initial data of the Einstein equations ([8])-([9]) 
satisfyinp the constraints (nnn-iiiD, then there exist unique classical solutions 
Pab and kab on [to, oo). 

In the above lemma the constant C/ is not assumed to be small. In the 
next section we will show that certain asymptotic behaviours are obtained by 
assuming smallness. 
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5.2 Bootstrap argument 

The argument which will lead us to our main conclusions is of a bootstrap type. 
One has a solution of an evolution equation and assumes that an appropriate 
norm of the solution decays in a certain rate. In our case the relevant function 
is F. Let with T < oo be the maximal interval on which the solution 

has the prescribed time decay rate. One improves the decay rate such that the 
assumption would lead to a contradiction, and concludes that the solution has 
the desired decay rate globally in time. 


5.2.1 Bootstrap assumption 


The bootstrap assumption is motivated by the fact that due to [11|, [1^ |18[ we 


expect the solutions to behave as solutions to the Einstein-Vlasov system. We 
will assume that the shear F satisfies on an interval [to,T) 


F < Cf(1 + t) 2, 

where Cf is a positive constant which will be determined later. 


(37) 


5.2.2 Estimate for the Hubble variable H 

Since we are assuming dSll), we have from (j36|) that 

S<CCft-V, (38) 

where S = g°‘^Sab = 'n°‘^Sab- Then, we use the estimate (1^ for F[ to get 

AnS 


3i72 


< CCft-^. 


(39) 


Integrating (I22L we have 




|t + / 3 l + 


where 


On the other hand, using the bootstrap assumption on F and the estimate (IMl) , 
we obtain 


/ < C{Cf + C/), 

where C is a positive constant depending only on the time origin to- Thus, we 
have 


TT > ±-L-l _ f _ 

“3 l + C{CF + Cf)t-^' 

which implies that we obtain a better estimate for FI, namely: 

H=\t-^{l + o({CF + Cf)t-^)). 


( 40 ) 
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5.2.3 Estimate of the shear F 


Our goal is to obtain an improved decay rate 

< C'F(l + t)”^+^ (41) 


on the interval [to^T) where F satisfies (l!?71l . More precisely, we will show that 
if F satisfies dSll), then there exists a small number 0 < <5 < | satisfying dm, 
where the small number 5 does not depend on T. This implies that F satisfies 
dUD globally in time, and completes the bootstrap argument. 

We now show that F satisfies (j^Tl) . If this is the case, there is nothing more 
to do. Let us suppose the opposite, namely that for any 0 < 5 < -g, there exist 
<1 and 0 such that to < ti < t 2 < T and 


Fit) <CFil + t)-^^\ t€[to,h], 

Fit)>CFil + t)-^^^, tG[ti,t2). 


(42) 


We consider the second time interval [ti,t 2 )- The evolution equation (IMl) is 
written as 

F _ - -F - — - ——^JF, 

and the quantities on the left side are estimated as follows: 

< jCf(l + t)-3 <CCf, 

M < CC,t-i < GG,. 

^<Y—- ^<CCft ii + t) ^<CCfC^ , 


where in the last inequality we have used (021) and the Cauchy-Schwarz inequal¬ 
ity. We now assume that the constants Cf and Cj are small. Note that in the 

last inequality although Cp ^ might be a big quantity, we have Cf as a factor as 

well which is independent and can be chosen such that the last quantity CfCp ^ 
is small enough. As a consequence we can find a small positive number (5i such 
that the following differential inequality holds: 


< {-2 + 6i)t-^F. 


Integrating the above differential inequality on \ti,t] with t G [ti,t 2 ), we obtain 

/ , \ -2+Si 

F(<)<C^(l+ti)-'+'(^-j 


which holds on [ti,t 2 )- Here, we have used E(ti) < (7^(1-I-ti)”^'*''^. Since 5 has 
been chosen arbitrarily and (5i does not depend on 5, ti, and t 2 , we may assume 
that S has been chosen as (5 = i5i. Hence, we obtain 

Fit) < CFil + t)-^+^\ 

which holds on [ti,t 2 ), and this contradicts to (1421) . Consequently, we have 
shown that there exists a small number 0 < (5 < 4 such that the shear F 

5 

satisfies (HIT) globally in time. 
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5.2.4 Optimal estimate of F 

We want to improve the estimate of F. For this reason we need an inequality 
in the other direction. From (1241) we have 

F > -3HF + 


Implementing now the estimates of F and / coming from the results of the 
bootstrap argument, we estimate the last term as 


< ^SFi 

\ + C{CF + Cf) 




-l+§ 


< CCfC^il + Cf + + t) 


-1+f 


< CCfC^t-^+^, 


where we used (gni), dMi), dm), and the fact that the constants Cf and Cf are 
small. We obtain 

F > -SHF - CC'/C|t-¥+i, 
and together with the upper bound H < 

F > -2t-^F - CCfC^t~V-+i, 


from which we obtain 

1 

t^F{t) > tgF(to) - CCfC} / s-i+Us, 

Jto 


and since ^ < -g, 

F{t) > (tlF{to)-CCfCiy-\ 


Here, F(to) should not be zero, and we take C/ small enough such that 

F{t) > hlF{to)t-^. (43) 

With this lower bound of F, we come back to the evolution equation (EH) and 
estimate it again as follows: 


< cCf{i +< CCFt-^^^, 

o 

8fS 9 


47r SabCf°-^ 

SH^ F 


< c4^F-i < CCfF{tQ)-h--s, 
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where we used (HTl) . (IMl) . and (H51) . We now obtain 


^ ^ -^‘~‘ i + c(cp + c,r. (‘ - + c, + C,r(,„)-J)ri) F 

< -2t-i(l - c{cf + Cf + 

where we used the first inequality of (HHl) , and consequently obtain an improved 
estimate: 

Fit) < CFito)t-^. 

Let us summarize the results obtained in the bootstrap argument part: 

Lemma 7. Suppose that ^ 0 is sufficiently small. There exists a small 

positive constant Cf such that if f is and satisfies (1351) . then there exists a 
small positive constant e such that the following estimates hold: 

H{t) = ^t-i(l + 0(et-i)), 

Fit) = 0{et-^), 

where the constant e depends on F{to) and Cf. 

Proof. If the initial data is small, then by continuity there must be an 

interval on which F satisfies Fit) < (7^(1 + t)~^ for some Cp- Let [toi?’) be 
the maximal interval, where T may be infinite. Suppose that T < oo. Since 
Fito) is small, the constant Cp can be chosen small, and we have shown in this 
section that in fact Fit) < C'f(1 +< C'f(1 + t)~i on [to,T). Now, by 
continuity F(t) remains less than Cpil+t)~^ for a short time after T, but this is 
a contradiction. Thus, T = oo, and the estimates obtained hold globally in time. 
The estimate of H and the improvement of F follow using these estimates. □ 


5.3 Estimate of the metric 

In this part we estimate the metric gat to obtain its asymptotic behaviour. In 
the previous lemma we only have asymptotic behaviours of the Hubble variable 
and the shear, but this is not enough to study the Boltzmann equation in the 
next section. To obtain global-in-time existence of solutions to the Boltzmann 
equation, we need an asymptotic behaviour of the metric at late times. We first 
estimate the metric gab and g°‘^, and then show that the scaled version g°‘^ is 
almost constant in the sense of quadratic forms. 


5.3.1 Estimate of the metric 


From the definitions m and (EHl) we have that 

9ab — 2 —t ^ gab ^(Tab- 

Using the usual matrix norm we obtain that 


||5a6(t)|| < ll5a6(to)|| +C 




(o-cdO-‘'‘^(s))2^||ga6(s)||ds, 
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and with Gronwall’s inequality we obtain 


||5at-(0ll < llffab(to)||exp H{s) - ^ + (CTcdCT°‘^(s)) = ds^ < C\\gab{to)\\, 


where we used Lemma [7] as follows: 

2 


H{s) - 


< Ces 

{<Jaba^\s))i = 2HF"^ < C^Es-^. 


Therefore, gab is bounded for all t > to, and we have 

\t~^gab\ < C, 

where we can choose the constant C as a constant which is independent of e. 
Note that HcTabU < Ugabll by [i^. Then, we can conclude that 

Ikabll < Cy/et~^, aab = 0{y/£t~i). 

Looking again at the derivative of gab and putting the facts which have been 
obtained together, we see that 


gab = ^). 

This is enough to conclude that there exists a constant matrix Qab defined by 

/ OO 2 \ 4 

2yH{s) - -s~^jgab{s) - 2s~'Saab{s)ds. 

Note that 

/ OO 

s~'^ds < Cy/et~^, 

and in this sense we obtain 

gab = (Sab + 0{^/£t~^)y 

We take the opportunity to note that in [l^ there is a mistake, where instead 
of 0{t~^) there is a wrong 0{t~^) in (61) and (62). By similar computations 
we can see that there exists a constant matrix such that the inverse is 
estimated as follows: 

We may choose a small number e again, if necessary, and combine the above 
results with Lemma [7] to obtain the following proposition: 

Proposition 1. Suppose that J^(to) ^ 0 is sufficiently small. There exists a 
small positive constant Cf such that if f is and satisfies ([55]). then there 
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exist a small positive constant e and constant matrices Gab and G°‘^ such that 
the following estimates hold: 

H{t) = |i-i(l + 0(£i-i)), 

F{t) = 

9ab{t) = t+i[gab + 0{et~^)^, 

g-\t) = + 

where the constant e depends on o,i^d Cf. 

5.3.2 Estimate of the scaled metric 

In this part we estimate the scaled metric g°'^ defined by g°'^ = t~^g°‘^. From 
algebra we know that in the sense of quadratic forms i.e., 

for any p* the following holds: 

cr°‘^PaPb < g°'^PaPb- (44) 

Introduce now for technical reasons £i and consider the following expression 
where we use the fact that = 2 /c“^ = 2 (ct“^ — Hg°‘’^): 

+ 1^-1 _ 2Hy-^'^g‘^^ + 2t-^^+^a‘^\ 

Using expression (HU) we obtain 

< (-eif-i + Cer^)t-^^r^PaPb + At-^^HF^g'^^paPb. 

Using the results of Proposition [U we have 

^{t~’^^g‘^’’paPb) < (Cy/et~'^ - £lt~^^{t~^^g°■^paPb)■ 

We can choose a positive number £i appropriately such that the following in¬ 
equality holds: 

(45) 

where the constant Si depends on the constant e of Proposition[I] which in turn 
depends on F{to) and Cf. Note that ei can be chosen small, since the constant 
e is small and t~^ < at late times. This estimate can also be stated as 
d{t~^^g°‘^)/dt < 0 in the sense of quadratic forms, and will be used to estimate 
the distribution function of the Boltzmann equation. 

5.3.3 Estimate of the determinant 

We apply (E51) to the results of Proposition [T] Integrating (1^ , we get 
det g(<) = det (/(to) exp ^6 f H{s)ds'j, 
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and applying it to Proposition [TJ 

det g{to)e~^^ < det g{t) < det g{to)e^^ . 

Since e is a small number, we may simply write 

^ detgjt) ty 

2 \to) - detg{to) “ \to J 

Applying this to the estimates 0 , we obtain 

9 ii 911922 — (912) 

From the explicit expressions of and in m and we derive 

\et\<Ct-i, \e^\<Cti, 
where the constants C depend only on initial data. 


( 46 ) 


(47) 


6 Existence of solutions of the Boltzmann equa¬ 
tion 


In this section we study the Boltzmann equation. We will assume that a space- 
time is given and will show that classical solutions to the Boltzmann equation 
exist globally in time. The spacetime will be assumed to satisfy the properties 
of Proposition[T]at late times. To be precise, we assume that the metric satisfies 
the following estimates. 


(A) Assumption on the spatial metric. There exists a positive constant 
C such that the metric gab and its inverse 5 “^ satisfy \gah\ < C and \g°'^\ < C 

4 7 4 L 

for gab = t^gab and g°'° = t~3g°'°. In an orthonormal frame, the metric can 
be written as gabec^\ = Vcd or gab = O^e^rjcd, and they satisfy |e)(| < Cfi 
and 10^1 < Ct3. Moreover, there exists a small number 0 < ei < such that 
d{t~^^g°‘^)/dt < 0 in the sense of quadratic forms. 


Under the assumption (A) on the metric, we will prove global-in-time ex¬ 
istence of solutions to the Boltzmann equation. The following inequalities are 
direct consequences of Lemma 01 


Lemma 8. Suppose that the metric satisfies (A). Then, the following estimates 
hold: 


\dp,p°\<Ct 3 , \dpM< 


CqH 


0+-4 


\dpiyfs\ < 


CqH 


,o+-4 


\dpM < c{p°)Hq^)h-i, \dp^^s\ < cip°)Hq°)h-i, 




< 


Cn^t- 


V'(n0)2 _ 


nae'^uj"nj 


ncW- 


< 


C 


-JnfffC' 


(48) 

(49) 

(50) 

(51) 


where the constants C do not depend on the metric. 
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Proof. The quantity in (15(11) can be written as 


d 




P° ^A+W 

and then the assumption (A) gives the first inequality. The other inequalities of 
are similarly proved. For the first inequality of (I49L we have 


<Ct ^\p-q\ 


p* g* 



pO qO 


Vl + |p|2 Vl + |g|2 


= Ct AiPa- <?a)(p“ - <?“) < Ct 3h^/Af, 
where we used (1281) in the last inequality, and then from (EH) we obtain 

^0 


\9pM = 


iy_ _T_ 

pO qO 


< C{p°)^ 3. 


The second estimate of (1491) is similarly obtained, since h < y/s, and the in¬ 
equalities EH and (IHT1) are also similarly proved from Lemma 01 □ 

Lemma 9. Consider a momentum p' in (na. The following estimate holds: 

\d,,pi\<Cp°{qy, 

where the constant C does not depend on the metric. 

Proof. Recall that p' is written in the representation (na as 


, Pi + qi h 

Pi = -^r- + - 


n^giae^i^ 


2 


and then dp^p'^ is given by 


n°g^aeA 


h (9p,p°)5iae^C*' 


d I Ah _ 

* 2 2 AAAAnAAY ^ V^(n0)2 _ 


h n°g^aeli'> 

2 (n°)2 — {ncC'^^'^Y 


dpu \j - (rice'^dA- 


We now estimate each term separately. The first term is clearly bounded by 
1/2. For the second term, we use (15^ (15^1) to obtain 


{dp^h)n^ g,aeA 




< cA ApAA Hp^ + g^)\9^aele\^A 




<c(p°{qA + Af), 

where we used (A). For the third term, we use (1571) . 

h {dpCAhAA 




and (1551) to get 

_2 ^/pA\g^aeA\V^ 


< cr^ ^ ivr < (J^O 
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For the fourth term, we use (EZD, (ESI), and EOD to get 


^ vW(/+l°)l5.ae?el9^/+ 9°)^ 


<C(g°) 2 (/ + gO+ 


0\2 


pu 

We combine the above estimates to obtain 

|5p,p'| < 

where we used the fact that > 1 and > 1 , and this completes the proof. □ 
Lemma 10. Consider a momentum p' in dni). The following estimate holds: 


ISp.p'l < C 1 


p\ p\ {P°? 
|u| + \h\ + |uP 

where the constant C does not depend on the metric. 

Proof. Recall that p' is written in (IT^ as 


(q 


0\3 


P^ = 


Pi + Qi 


+ n ( 9iae-luj + ( — 1 


Ugeluj'^ni 

g^^'^ncUd 


and directly differentiate this with respect to pk to obtain 


a„p: = I + 


giae-b^ + ( — 1 


ngclw^ni 

g^'^UcUd 






^-1 


Uge'^uj'^ni 


g’^’^UcUd 


Note that in an orthonormal frame, 




Pk 


Uae'^uj^rii 


_ g'^'^HgUd 


nae'^LO''n^ 


g^'^UcUd 


\hbuj^e'ti 


< Ct^. 


We now estimate the above quantities separately. The first term is bounded by 
1 / 2 , and for the second and the third terms we use (HRl) and ( 1 ^ to get 


{^Pk a b 

--gzaCfeW 


< 


Cq° 


and 


jdpkh) f n° _ \ ugefuj' 

2 Va/s J g^^nc‘ 


fh 

nd 




UgC'^uj'^ni 


hy/s g’^’^UcUd 


Cq^n^ 

~1p 


< 


For the fourth term, we note that 

nO 


d, 


'Pk 


^-1 




° n° 

~^Pk V^j 


dpkP' ^ n 
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and then we can use (H51) . Then, we have 


-a. 


Pk 


^-1 




nae^uj'^ni 


g’^^ncUd 


< Ch 




For the fifth term, we use (EZl) SjUcI (151 1) to §Gt 




ngegcu^n, 

_ g^’^ncUd _ 


< 


Cn° 


y/n^ 

We combine the above estimates and apply n° < and 

\d,,p[\ < C(^l + ^ ^ (g°)3, 

and this completes the proof. 


to get 


□ 


Now we are in position to construct classical solutions to the Boltzmann 
equation. We choose the weight function as 


exp(< ’^^g°‘^PaPb) = exp(f 


“Ei + 7? |.A|2 


i^n, 


and define 


\\f{t)\\ = ^sup ||exp(T ^^g°‘^PaPb)dpJ{T,p^)\ : 0 < r < f, p* e R^j, 

n=b,i 

where ei is the constant given in the assumption (A). 

Lemma 11. Suppose that a spacetime is given and the metric satisfies the 
assumption (A). If f is a solution of the Boltzmann equation, then f satisfies 
the following estimate: 

exp(f“®ig“Vp&)/(^,P*) < \\fito)\\ +C'||/(t)f. 

Proof. To the equation in the form of m we multiply the weight function to 
obtain 


dt 


exp(f ‘'^g°‘^PaPb)f{t,p*) 

d 


= ^^g‘'^PaPb)exp{t ‘^^g‘'‘^pcPd)f{t,p.,) + exp{t g‘'^PaPb)Q{f, f){t,p.,) 

< {detg)~i JJ VMcr{p,.,q,.,U})f{p'Jf{qi)exp{t~^^g°^^PaPb)dujdq,,, 


where we ignored the loss term Q- and used the assumption (A) such that 
d{t~^^g‘^^PaPb)/dt < 0. We now use LemmalHlas follows: 

f{p*)f{Q'*) exp{t-^^g°-'^paPb) 

= (exp(f-^ig“*'p>g)/(p'j) (^exp(f-®^5“*'g;gg)/((7')) 

X exp(f"®ip“'’(pgp6 + qaqb - p'aP'b - <lWb)) exp{-t~'^^g°-^qaqb) 

< c(^exp{t--^r%Pb)fip:)){^Mt-^^r\Wb)fiQ'*)) expi-t-^^r\aqb) 

< C\\fit)f exp{-t~'^^g°-^qaqb)- 
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Since vm and (t(p*, g*,a;) are bounded, we have 

(det5)“5 JJ VMcr{p^,qt.,u})fip'Jf{qi)exp{t~‘^'^g‘^’’paPb)du}dq^ 

< C{detg)-^\\f{t)f J exp{-t-^^g°‘^qaqb)dq^ 

= C\\f{t)f J \qf)dq < C\\f{t)ft- 


|2 .- 2 +^£i 


Since 0 < £i < is small, the quantity t ^+ 2^1 is integrable, and we obtain 


that: 


exp(t ^^ 5 “ PaPb)f{t,P*) < ||/(io)|| + C\\f{t)\Y 


and this completes the proof. 


□ 


Lemma 12. Suppose that a spacetime is given and the metric satisfies the 
assumption (A). If f is a solution of the Boltzmann equation, then f satisfies 
the following estimate: 

|exp(t"^ig“Wb)9p,Q(/,/)| < C'||/(t)ft“2+5^b 

Proof. In the proof of this lemma we need to consider the second expression of 
the Boltzmann equation (HU- We first consider m and take derivative to it 
with respect to pk- Multiplying the weight function, we obtain the following 
quantities: 

exp(t“®^g“Vp&)5p;,Q(/, /) 

= {detg)~i JJ 9^, g*, w) (^f{p'.,)f{qi)-f{p*)f{q*)^exp{t~^^g°'^PaPb)dujdq., 

+ {detg)~^ JJ VMcr{p*,q*,u})dp^ f{p'*)f{q*) exp(t“^ip“VaPb)'^wdg* 

- {detg)~^ // VMcrip*,q*,u})dp^f{p.,)f{q.,)exp{t~^^g°‘’’paPb)du;dq^ =: Ji + J 2 + Ja- 


The estimate of J 3 is easily obtained and is 

\J3\<C\\fit)\\H-^+i^K 

For Ji, we note that the scattering kernel is bounded and 

< ^\dp,h\ + < Cq°t-i < Cq°, 

where (HZl) and has been used for the first and the second quantities, and 
dUl) and (H51) for the third quantity. Following the same calculation as in Lemma 
m we can see that Ji is estimated as 


|Ji| < C(detg) "^\\f{t)f J q° exp{-t g°‘'’qaqb)dq^ 

= c\\f{t)r J y/rTi^exp(-t--^+t|9ndg 


= c\\fit)rt 


2 I —2+4ei 


1 +f-3 


-o+El 


z| 2 exp(-|z| )dz < C||/(t)|| t 


2 .- 2 +^ei 
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where C may depend on initial data. 

We now estimate J2- Since the derivative term is expanded as 




the quantity J2 is separated into two integrals. The calculations are exactly 
same, and we only consider the first case, i.e., J2 = J21 + J22, where 


J2i = (detg) JJ VMfT{pt:,q^,u})idpJ){p'J{dp^Pi)f{qi)exp{t g°-^PaPb)du;dq^. 


Note that the integral J21 can be written as 

J 2 i = (detg )“5 / VMKexp(t~'^^g‘^^paPb)dq^, 
JR 3 


where 

K= f a{p^,q^,uj)(dpJ){p'J{dp^p'i)f{q'Jduj. (52) 

The second form of the Boltzmann equation m is applied to the quantity K. 
By the transformation cc 1—^ with the Jacobian (1591) . we obtain 


K = 


cr{p*,q*,0 {dpi /) (p* ){dp^p'i)fiq'*) 


y/s{n°)‘^d^ 


{{noy - 


(53) 


where and q'^ are now understood as m- The estimates of J21 are separated 
into three cases as follows: (a) |pp < 1, (b) |pp > 1 and |pp < djgp, (c) |pp > 1 
and IpP > 4 |gp. 


(a) In the first case, p° is bounded, and we can use dMI) with the estimate of 
Lemma O 

P° = Vl+W < V2, \dp,pi\ < Cp°{qy < C{qy. 

Together with (l?71) and the integral K in (1551) is estimated as follows: 

\K\<C [ (7{p^,q^,0\dpJ{p'J\{q°ff{qi)dC 

Js^ 

The integral J21 is now estimated by the same argument as in the previous 
lemma: 

\j2i\<c\\fmH-^+^^^K 

(b) In the second case, we use again the second form of K in (l53ll . Then, we 
have 


p° = < ^T + W < 2g°, \dp,p',\ < Cp\q^f < C{q^f. 

Together with (1571) and (1551) . we have 

\K\<C ( (T{p^,q^,i)\dpJ{p'^)\{q°yf{q'^)d£,, 

and therefore again 

\J2i\<C\\f{t)\\h-^+^^^K 
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(c) In the third case, we use ( 15 ^ with the estimate of Lemma [TUI Note that 


/ = \/l + \P? < \/2|p|, \p±q\> \p\-\q\>^\p\. 


Hence, Lemma [TUI implies that 

< C{q^f. 

The quantity K in (IK^ is estimated as 


\K\<C 



cr{p*,q*,u})\dpJ{p'J\{q°ff{q'Jduj 


and by the similar calculations we obtain 




We combine the above estimates to get 


|exp(t ^^g’^^PaPb)dp^Qif,f)\<C\\f{t)\\h ^+2^1 


which completes the proof. 


□ 


Lemma 13 . Suppose that a spacetime is given and the metric satisfies the 
assumption (A). If f is a solution of the Boltzmann equation, then f satisfies 
the following estimate: 


|exp(t ^^g°‘^PaPb)dpJ{t,p^)\ < ||/(io)||+ (^11/(011^■ 


Proof. For simplicity, let us write 


w = t ’^^g'^'^PaPb, (l) = e'^dpj, M = e'^dp^^QifJ). 

Then, the Boltzmann equation can be rewritten as 

dtfi = {dtw) 4 > + M. 

Note that w is nonnegative and decreasing by ( 1 ^ . Multiplying e to the 
above ODE and integrating it from to to t, we have 



and then 



Applying Lemma [T^l we conclude 

\exp{t-^^g‘^^PaPb)dpJ{t,p.,)\ < \\f{to)\\+C\\f{t)\\^, 


which completes the proof. 


□ 
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It is now easy to prove the global-in-time existence of classical solutions to 
the Boltzmann equation. We combine Lemma [TT] and [T 3 ] to obtain 




and by the well-known arguments as in I, i i [I^ fill 1^ . we obtain the 
global-in-time existence of classical solutions for small initial data. 


Proposition 2. Suppose that a spacetime is given and the metric satisfies the 
assumption (A). There exists a small positive constant Eq such that if f{tf) 
is and satisfies ||/(to)|| < ^O; then there exists a unique classical solution 
to the Boltzmann equation corresponding to the initial data /(to) satisfying 
||/(t)|| < C'll/(to)|| on [to,oo). The solution f is nonnegative and satisfies in an 
orthonormal frame 

f{t,p) < C'||/(to)|| exp(-t"''i+5|p|2), 


where Ei is the constant given in the assumption (A). 


7 Main results 


We are now ready to prove global-in-time existence of classical solutions to 
the Einstein-Boltzmann system ©-(ini). The existence is proved by a stan¬ 
dard iteration method. Suppose that initial data gab{to), kabftfi), and /(to) are 
given such that P(to) 7^ 0 is sufficiently small, and define an iteration for {gn}-, 
{k„}, and {/„} as follows. Let (go)ab(t) = tagabih) and (fco)ab(t) = kabih). 
Then, there exists (6»o)b satisfying {go)ab = {do)l{ 0 o)'^r]cd, which is given by 
( 0 o)b (t) = t^ 0 ^(to), and let (60)5 be the inverse of ( 0 o)b • Here, ^5 (to) is uniquely 
determined by (5o)ob(to) = (^o)a(to)(^o)b(to)»?cd- Note that (go)ab(t) = 5ab(to), 
be., {go)ab is a constant matrix, hence clearly we have d{t~'^^{go)°'^)/dt < 0. 
In other words, go satisfies the assumption (A) of Section [51 Then, by Propo¬ 
sition [21 there exists a small positive constant eo such that if ||/(to)|| < £0, 
then there exists a unique classical solution /o, which is the solution of the 
Boltzmann equation in a given spacetime with metric go and satisfies ||/o(i)|| < 
C||/(io)||- Now, suppose that /„ is given such that ||/n(0ll — C'||/(<o)|| with 
II /(to) II < £o- Note that, since Ei < ^, we have in an orthonormal frame, 
fn{t,p) < C\\f{to)\\exp{—ti\\p\'^). We now apply Proposition |T] Since we have 
assumed that ||F(to)|| 7^ 0 is sufficiently small, there exists a small positive 
constant C/ such that if /„ satisfies (jSS]), then gn+i and exist globally in 
time, which are the solutions of ODEs, which result when g and k of di)-® 
are replaced by gn+i and fc„+i, respectively, and p and Sab are constructed by 
fn via (fTBl) and (EZD, respectively. Taking £0 smaller, if necessary, we can con¬ 
clude that if fn is given such that ||/ra(t)|| < (711/(^0)11 with ||/(to)|| < £0) then 
gn+i and kn+i are constructed such that they satisfy the assumption (A). Now, 
Proposition 121 again applies, we obtain fn+i, and this completes the iteration. 

We have obtained iteration functions {gn}, {kn}, and {fn}- The estimates 
of Section [ 5 ] show that the following quantities are uniformly bounded: 


t^Hnit)--t, \t^Fnit)\, \ign)abit)\, | (5„)“''(0 1, 1^" (ffn)ab (i) |, \t\kr{t)\, 
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where Hn and Fn are understood as the Hubble variable and the shear con¬ 
structed by gn and kn- Using the evolution equations ([S])-®, we can see that 
the following quantities are also bounded: 

\t~Hkn)abit)\, \tiikn)abit)\. 

The estimates of Section [6] show that the following expressions are bounded: 

exp{t~^^{gnT''paPb)fn{t,p*), | exp(t“®i/„(t,p*)|. 

Taking limit, up to a subsequence, we find continuous functions g, k, and /. Us¬ 
ing the equations dSl)-® again, it turns out that g and k are classical solutions 
of ([Hl)^®- Similarly, we can show that / is a classical solution of the Boltz¬ 
mann equation with additional estimates. Uniqueness is also easily proved by a 
standard argument, and asymptotic behaviours are obtained by the arguments 
of Section [S] and [51 The following is the main theorem. 

Theorem 1. Consider the Einstein-Boltzmann system with Bianchi I symme¬ 
try @-( 1171 ). Suppose that the assumption on the scattering kernel holds and 
the Hubble variable is initially positive, and set the time origin as in (|TO)) . Let 
Pabito), kab{to)i OL^d /(to) be initial data of the Einstein-Boltzmann system sat¬ 
isfying the constraints (nnii-diii) such that /(to) is and U(to) 0 is suf¬ 
ficiently small. Then, there exist small positive constants Eq and Ei such that 
if II/(to)II < Eo, then there exist unique classical solutions gab, kab, and f to 
the Einstein-Boltzmann system corresponding to the initial data. The solutions 
exist globally in time, and the distribution function f is nonnegative. Moreover, 
there exist a small positive constant e and constant matrices Qab o,nd Q°‘^ such 
that the following estimates hold: 

H{t) = ^t-^{l + 0 {Et-^)), 

F{t) = 0 (£t- 2 ), 
gab{t) = t+^ (Oab + 0{Et~^)y 

g-\t) = t- 5 (a“'' + 0(et-i)). 

The distribution function f is bounded as ||/(t)|| < C'||/(to)||, and in particular 
it satisfies 

f{t,p) < C'||/(to)|| exp(-t“®i+^|p|2) 
in an orthonormal frame. 

With the results of the theorem, we are able to show that the spacetime 
behaves as dust at late times. We show that — tends to zero which means 

p , 

that the “pressure” tends to zero in comparison with the energy density. In an 
orthonormal frame the estimates of S and p are as follows: 


P= f{t,p){l + \p\r^dp<C exp(-t-"i+3|p|^)(i + |p|^) 2 dp<c< 


i-2-l-f El 


/R 3 


/R3 


\Sab\< [ f{t,p)\pa\\pb\{l + \p\^)-^dp<C [ exp{-t---+i\p\^)\p\^dp<Ct-^+i^K 
JR3 JR3 


31 


Since S = rj^-^Sab, we obtain that S = 0 {t 3“+ 2=1) and p = 0 {t ^+2^1). For 
—, we need an estimate of p in the other direction. We can use the fact that 
p > 7 V°. Using the upper bound of H in ( 1 ^ and the differential equation ( 1 ^ . 
we have 

iV° = > - 2 t-^N°, 

which can be integrated to obtain 

We obtain the following corollary. 

Corollary 1. Consider the solutions of the Einstein-Boltzmann system given 
in Theorem [IJ There exists a small positive constant e such that the matter 
terms have the following asymptotic behaviour at late times: 

p = 0(t-2+-), 

S' = 

- = 0 {t-i+^). 

p 

Another consequence of the theorem is the following: 

Corollary 2. Consider the solutions of the Einstein-Boltzmann system given 
in Theorem]^ The spacetime in consideration is future geodesically complete. 

In order to prove this, the main point is to consider the relation between 
the time variable t and the affine parameter r along any future directed causal 
geodesic which is: 


dt p° s/l+WpiFj 

The last inequality is obtained by realizing that the pi are constant along the 
characteristics for Bianchi I spacetimes and that the inverse of the metric decays. 
Integrating the inequality we see that r goes to infinity as t does. 


8 Outlook 


We have shown under small data assumptions and an assumption on the scat¬ 
tering kernel, that solutions to the Einstein-Boltzmann system with Bianchi 
I symmetry are future geodesically complete, that they isotropize and have a 
dust-like behaviour at late times. Spacetimes of Bianchi type I are the sim¬ 
plest case of the Bianchi class A. It would be thus of interest to consider other 
Bianchi types. Here we have treated massive particles and a natural general¬ 
ization would be to consider massless particles or to consider a collection of 
particles with different masses. Our assumption on the scattering kernel is a 
technical assumption and we plan to remove or at least improve this assumption 
in the future. In the present paper we were dealing in a homogeneous setting. It 
is a simplifying assumption which should be removed as well, since conclusions 
do not have to prevail in an inhomogeneous setting as showed in 231 . 
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Appendix 

In this part we study the Boltzmann collision operator and parametrization of 
post-collision momenta in relativistic cases. We first consider the Minkowski 
case, in which case two different representations of the collision operator are 
known, which are given by Glassey and Strauss Q and Strain [^. We use the 
argument given in Q to derive a similar form of the collision operator, and it 
will be shown that this form of the collision operator can also be derived in a 
direct way. To extend these representations to the Bianchi type I case we apply 
the orthonormal frame approach. An explicit expression of orthonormal frame 
will be given, and then the representations for the Bianchi type I case will be 
obtained. 

Collision operator in the Minkowski case 

Let us consider the collision operator of the Boltzmann equation in the special 
relativistic case. In its original form we have 

^ W{p,q\p',q^)(^fip')fiq') - /(p)/(g)) 

where the quantity W is called the transition rate [ 2 ^, and it is given by 

W (p, q\p', q') = ^a{h, 9)6^^^ {p^ + q^ - p'^ - q'^), 

where is the four-dimensional Dirac delta function which expresses the con¬ 
servation of energy and momentum. It is known in that the nine dimensional 
integration above is reduced to 

and the post-collision momentum is parametrized as follows: 

/ p^ + q° h {n ■ uj) \ 

2 ^2 yi 

V 2 2^^ 

where w G To obtain the above representations, a suitable Lorentzian 
transformation should be considered. 

On the other hand, the nine dimensional integration can be reduced in a 
different way without using Lorentz transformations so that a different form of 
the collision operator is derived. The post-collision momentum may be assumed 
to have the form of p' = p -|- as in [ 7 [ , but we replace this by 

= C e ( 56 ) 

and then follow the calculations of Q . We first consider the integration over K.® 


UJ 


n 


- 1 


(n • a;)n* 
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with respect to p' and q': 


IIrb + q^-p''" - q'^)(^f{p')f{q') - fip).fiq)) 

= a{h, 0)6{p° + g° - - g'O) (f{p')f{q') - f{p)f{q)) ^ 

= 0)5{p^ + <z° - p'° - q'^) [f{p')f{q') - f{p)f{q)y-^, 

where q’ = p + q — p' in the first identity, and we assume (1551) in the second 
identity. For the Dirac delta function we use the property S(x) = |a|(5(aa;) to 
obtain 


d{p° + q°- p'° - q'°) 

= (p° + + q'°)s(^ip° + qy - (p'° + g'°) 2 ) 

= 2(p0 + g>((pO + - (P'° + O") • 

Similar arguments give the following calculations: 

2{p° + q°)d[ip'°+q'y-ip°+q°y 

= 2(p0 + q°)6(2p'\'° - (p° + qy + {p'y + (g'°)2) 

= 8(p° + q^)p'\'H[^{p'y{q'y - ((p° + qy - {p'y - , 

and let p = p{r) denote the quantity contained in the Dirac delta function. 
Then, it is rewritten as follows: 

pir) = Ayyyy - y + qy + 2y + qyyy + yy) - yy + y)y 

= -y+qy - yy - yyf+^y+qyyy + yy. 


Note that 

yy-{q'yr = i\pf-\q'\y = 


p + q 


+ 


P + q 


- 


‘2^2 


= (2(p + g) -y =4(n-0 r 


and 


{p'y + {q'y = 2+ 


p + q 


+ 


P + q 


- 


= 2 + 


2 


+ 2r^ 


where n = p + q and = p^ + q^. Therefore, the quantity p{r) is written as 

p{r) = -yy - (yy - y)y + 2{nyyy + {q'y) 

= -{ny - 4(n • yr'^ + 2{ny (2 + + 2r^^ 

= r^{4{ny - 4(n • 0") - {ny{{ny - |np - 4). 


34 
















We may write (n°)^ — (n • for t°' = {n ■ n°^) to obtain 

p{r)=A{t^ey -h^iny, 


where we used (n°)^ — |np — 4 = s — 4 = h?. We plug this to the above 
calculation: 

lo X “ P'° ~ 

= 8^ X 9°y°9'°'^(p(^))(/b')/('?') - f{p)f{q)) ^pfo^^Q 

= i I a{h,9)n°6(4{tc,t°‘)r‘^ - (^fip')f{q') - f {p) f d^dr 

L ^-/(»>)/(«))’■“■«* 

' ®i i. ■ 2^) (r(i'')/(<,') - mmy^Jidr 

‘L L “2^)- fMMydfdr 

=U. 

Since tat°' = (n°)^ — (n ■ we combine the other quantities to get 


Q(fJ) 


hs(n^ya(h, 9) 


/s2 4p°(70((n°)2 — (n • ^)2)3/2 


{f{p')f{q')-f{p)f{q))didq. (57) 


h ® 

From the above calculations we find r = to write the post-collision 

momentum as 


/ p° + q° ^h (n-g) \ 


o 

1 

2 2y'(n0)2 -(n.O2 

1 

-K g'= h rPSy^ 


V 2 2 y^(n0)2 - (n • 0" / 


(58) 


where ^ G S^. 

We may now compare ([Si]) and (EB to see that the two parameters w and ^ 
have the following relation: 


duj 


\/i(n '’)2 

((„0)2 _ („.^)2)3/2 


d^, 


(59) 


and this can be proved by direct calculations. Suppose that a post-collision 
momentum p'“ is given. It has two different representations (1551) and (1551) . and 
we have the following identities: 


n ■ oj n ■ y, 

yr = V^p)^ (n • ’ 

nP \ (n • u})n^ 

7i " V |nP " v'(^0)2-(n.#' 
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Taking coordinates satisfying n = (0, 0, |n|), we write oj = uj{a,j3) and ^ = 
9) as follows: 


w = (sin a cos/3, sin a sin /3, cos a), 
^ = (sin^cos6*,sin(/)sind,cos(/()), 


and then the second identity reduces to 

/ . ^ \ n°(sin(()cos0,sine/)sind,cos(()) 

sin a cos p, sin a sin p, —cos a = - , =-, 

V - jnl^ COS^ (/> 

where «,</>€ [0,7r] and 13,9 G [0, 27r]. This quantity describes a certain two- 
dimensional surface in and the above identity shows that the surface is 
parametrized in two different ways. This implies the surface element can be 
computed in two different ways. Let X denote the surface. Then, we have 


dS = \daX X dpX\dadl3 = \d(j)X x deX\d(l>d9. 
By direct calculations, we have 


( T^ \ 

COS a cos (3, cos a sin /3, — sin aj , 
df}X = ^ — sin Q! sin (3, sin a cos (3, 0^ , 


and 


( Tv Tv \ 

—= sin^ a cos (3, —= sin^ a sin /3, sin a cos a ). 
VS VS ) 

Hence, we have 

,2 _ 


lev -\.r ev tx" 12 ) • 4 -2 2 SID. Q^//n\2*2 2\ 

\docX X d^X\ = -sm a + sm a cos a =-(^(^ ) sin a + scos aj 

({ny - \y cos^ a) = ^ ((n°)2 - (n • c.)^), 


sin^ a 


and therefore we obtain 


|a„x X a,x\iaifi = 

VS 


(60) 


where we used dw = su\adadf3. We now consider the second case X = X{<p, 9). 
By direct calculations again, we have 


dd,x = a 


(sin (j) cos 9, sin (j) sin 9, cos </>) 


a/( n^)^ — |np cos^ 

H- , z (cos 6 cos 9, cos 0 sin 0, — sin 6) 

a/( n‘’)2 — |np cos^ (j) 

-n^\ysm(j)COS(j) a ■ ± ■ a 

- - —^(sin a cos 6», sm </> sm 6», cos </>) 

H- , _ _ , _ ^ z (cos (j) cos 9, cos a sin 0, — sin 4 >), 


a/ {ny — |np cos^ cj) 
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and 


dgX = =(—sin(/)sin0,sin0cos6>,O). 

y [vPy — |np cos^ (j) 

Note that x = \d^X\^\d 0 X\^, where the quantities are calculated as 

|a 2 _ (^°)^l^l'‘sin^<?^cos2(;() {nPf 

' ' “ ((n0)2 - |n|2cos2<;i)3 + (n0)2 - |n|2cos2</)’ 


and 


\deX\^ 


(n°)^ sin^ (p 
(n °)2 — |n |2 cos 2 (j)' 


Hence, we have 


la^X X dgX\'^ 

(nP)"^ sin^ (p 
((n °)2 — |n |2 cos 2 p)'^ 
(nP)^ sin^ p 
((n0)2 — |n|2 cos2 p)'^ 

(n°)^ sin^ p 


{{n°P - {n ■ 
{rP)‘^ sin^ p 
{{rPy - {n ■ 


^|n|'^ sin^ pcos^ p + {{rPY ~ cos^ pY^ 
^|n|^ cos^ p + — 2(n°)^|np cos^ p'^ 

(|np(n.0^ + (n'’)4-2(n°)2(n.ef) 


^ {n^YsiiPp / 0 2 _ s{n-0‘^ \ 

{{n°Y - {n ■ YYY y ’ (nO)2 -(n-e)V' 

On the other hand, we can use the hrst relation between ^ and uj to the last 
quantity as follows: 


yr- 


s(nyY 


(n 0)2 - (n • iP 


= K)"-(n-w )2 


Then, we obtain 


Icl^X X dgX\dpd9 = 


(713)2 y^(n0)2 - 
((770)2 -(n-02)3/2 


(61) 


where d^ = sin pdpdO. Finally, we know from calculus that the surface element 
does not depend on the way in which the surface is parametrized, and therefore 
we compare (pOl) and ((HTIl to get the desired result: 


dio 


y~s{rPY 

((770)2 -(71-02)3/2 


de, 


which we have already observed in (I59F 
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